Quantum Heisenberg antiferromagnets: 
a survey of the activity in Firenze 
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Over the years the research group in Firenze has produced a number of theoretical results con- 
cerning the statistical mechanics of quantum antiferromagnetic models, which range from the the- 
ory of two-magnon Raman scattering to the characterization of the phase transitions in quantum 
low-dimensional antiferromagnetic models. Our research activity was steadily aimed to the under- 
standing of experimental observations. 



I. INTRODUCTION 



The Heisenberg model may well be considered the cor- 
nerstone of the modern theory of magnetic systems; the 
reason for such an important role is the simple struc- 
ture of the Hamiltonian, whose symmetries underlie its 
peculiar features. The basic forces which determine the 
alignment of the spins, are represented by the exchange 
integrals J's. At variance with the ferromagnet where 
the parallel alignment is promoted, in the antiferromag- 
net a lot of peculiar arrangements of the spins can occur, 
with strong differences between classical and quantum 
systems. As matter of fact, also for nearest-neighbor anti- 
ferromagnetic interactions the ground state of the Hamil- 
tonian is different from the Neel state with antialigned 
spins, and the (staggered) magnetization shows the so 
called spin reduction with respect to the saturation value 
also at T — 0. The linear excitations of an antiferromag- 
net can be roughly associated in two families and pair 
excitations with vanishing total magnetization are possi- 
ble: the fact that the total momentum of these can be 
close to zero allows for their investigation by light scat- 
tering. 

While these peculiar features of antiferromagnetism al- 
ready occur in three-dimensional (3D) compounds, they 
are more pronounced in the low-dimensional ones, where 
other effects caused by the enhanced role of classical and 
quantum fluctuations are present and exotic spin config- 
urations associated with field theory models can appear. 
Indeed, the last two decades have seen a renewed inter- 
est both in the case of the one-dimensional (ID) quan- 
tum Heisenberg antiferromagnet (QHAF), for which a 
peculiar behavior of the ground state vs spin value was 
predicted 1 -, and of the two-dimensional (2D) QHAF, be- 
cause of its theoretically challenging properties and of the 



fact that it models the magnetic behavior of the parent 
compounds of some high-T c superconductors 2 ^. The ex- 
perimental activity on 2D antiferromagnets stems from 
the existence of several real compounds whose crystal 
structure is such that the magnetic ions form parallel 
planes and interact strongly only if belonging to the same 
plane. As a consequence of such structure, their magnetic 
behavior is indeed 2D down to those low temperatures 
where the weak interplane interaction becomes relevant, 
driving the system towards a 3D ordered phase. 

In addition, the 2D Heisenberg model can be enriched 
trough symmetry-breaking terms - we considered easy- 
axis (EA) and easy-plane (EP) anisotropy, as well as 
an external uniform magnetic field - which are useful 
to reproduce the experimental behavior of many layered 
compounds. In the EA case one is left with a discrete 
reflection symmetry and the system undergoes an Ising- 
like phase transition. In the EP case or when a mag- 
netic field is applied the residual 0(2) symmetry prevents 
finite-temperature ordering^, but vortex excitations are 
possible and determine a Berezinskii-Kosterlitz-Thouless 
(BKT) transition between a paramagnetic and a quasi- 
ordered phase. In spite of the tiny anisotropies of real 
systems (usually < 0.01 J), it can be shown that they 
dramatically change the behavior of the spin array al- 
ready at temperatures of the order of J. 

In this paper we report about the progresses in the 
theory of Heisenberg antiferromagnets that have been 
obtained by our group in Firenze. The early work on 
the theory of two-magnon Raman scattering is summa- 
rized in Section [H] while the following Sections report 
about the recent activity on low-dimensional antiferro- 
magnetism. Section IHII is devoted to ID models, and 
concerns the study of the effect of soliton-like excitations 
in the compound TMMC, as well as the anisotropic spin- 
1 model, for which a reduced description of the ground 
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state allows one to investigate the quantum phase tran- 
sition in a unitarily transformed representation and to 
obtain quantitative results for the phase diagram. Sec- 
tion llVl concerns the theory of the isotropic 2D QHAF, for 
which we reproduced the experimental correlation length 
by means of a semiclassical approach, also deriving the 
connection with (and the limitations of) famous quantum 
field theory results. In Section we summarize several 
recent results concerning the anisotropic 2D QHAF, with 
emphasis onto the different phase diagrams and the ex- 
perimentally measurable signatures of XY or Ising behav- 
ior. Eventually, in Section lVll results on the 2D frustrated 
J1-J2 isotropic model are described. 



II. TWO-MAGNON RAMAN SCATTERING IN 
HEISENBERG ANTIFERROMAGNETS 

The scattering of radiation is a very powerful tool 
to study elementary excitations in Condensed Matter 
Physics. Any complete experiment gives rise to a quasi- 
elastic component due to non-propagating or diffusive 
modes and to symmetrically shifted spectra correspond- 
ing to the states of the system under investigation with 
an amplitude ratio governed by the detailed balance 
principle. The most sensitive probes for this investi- 
gation are undoubtedly thermal neutrons, because the 
characteristic energies and wavevectors fit very well with 
those of the magnetic elementary excitations. However, 
light-scattering experiments can require a simpler appa- 
ratus and offer a better accuracy, although the transfer 
wavevector k is much smaller than the size of the Bril- 
louin zone so that usually only the center of this zone 
can be directly probed. In spite of this, two-spin Ra- 
man scattering involving the creation and destruction of 
a pair of elementary excitations can be performed, with 
the contribution of two magnons having equal frequencies 
and opposite wavevectors. This two-magnon scattering 
is expected to be spread over a band of frequencies in an- 
tiferromagnets. However, the density of states strongly 
enhances the contribution of zone-boundary (ZB) excita- 
tion:-: 5 , i.e., at fc~fc ZB . 

While in ferromagnets the two-spin process is only 
due to a second order mechanism, orders of magnitude 
smaller than the first order one, in antiferromagnets a dif- 
ferent independent process is permitted, stronger than 
the corresponding one for single-spin spectral. Specifi- 
cally, an exchange mechanism does not change the total 
^-component of the spins: exciting two magnons in the 
two different sublattices (AM = 0)£ is the dominant scat- 
tering process. 

The one-spin Raman scattering peak disappears at the 
Nccl temperature because it probes the smallest wavevec- 
tors, related with the long-range correlations. In con- 
trast, two-magnon Raman scattering essentially probes 
the highest wavevectors, related to short-range correla- 
tions. Therefore two-spin Raman scattering features per- 
sist also in the paramagnetic phased where short-range 



order is still present. 

Let us consider the following antiferromagnetic Hamil- 
tonian with exchange integral J > 0, z nearest neighbors 
with displacements labeled by d, and two (a, b) sublat- 
ticea&: 

T~L = 77 Si,g-Si+d,b ■ (1) 
id 

The scattering cross section S(to) turns out^ to be pro- 
portional to the Fourier transform of (M(O)M(t)), where 

M = J2 Mk S kS-k, (2) 

k 

is the effective Raman scattering operator. 

Many antiferromagnetic compounds can be mapped 
onto this model, even though a small next-nearest- 
neighbor exchange interaction without competitive ef- 
fects, as well as anisotropy terms could be present. For in- 
stance, there are 3D perovskite and rutile structures (e.g., 
KNiF 4 , NiF 2 ) and 2D layered structures (e.g., K 2 NiF 4 , 
LaCu0 2 ). 

Let us remember that the exact ground state is not ex- 
actly known, except in ID models with S = 1/2 or S = 00 
(i.e., the classical case): in the latter case it coincides 
with the 'Neel state' with antialigned sublattices. 

In the ordered phase the theory can be developed in 
terms of two families of magnon operators (a k , j3 k ), 
through the Dyson-Maleev spin-boson transformation 
and a Bogoliubov transformation: 

H = E + Ho + V , (3) 

where E is the ground state energy in interacting spin- 
wave approximation and 

n = J2^(a{a k +p{p k ) (4) 

k 

is the quadratic part of the Hamiltonian of a magnon 
gas whose frequencies, renormalized by zero-T quantum 
fluctuations, are 

u> k = JSz(l + £)y/l^rt; 7 fe = ~£e- ifed . (5) 

d 

The last term in the Hamiltonian, V, represents the 
four-magnon interaction, whose most significant term 
refers to two magnons of each family and turns out to 
be: 

V = 2 Jj~ S i+p,q'+p' I tf,p P > a q a q'PpPp' > ( 6 ) 

qq'pp' 

where the coefficients 1™^, , are known functions of 7*.. 

In the Hartree-Fock approximation^ the temperature 
dependent Raman scattering operator J2J can be written 

M = a(T) *fc(a*/3fc + 44) , (?) 

k 
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FIG. 1: Theoretical two-magnon spectra in KNiF3 at differ- 
ent temperatures--. 



with u)k{T) — a(T) cuk- The two magnons created or de- 
stroyed by the operator J7J interact through V as given 
by (JBJ, so that the peak of the cross section S(uj) appears 
at values smaller that 2w ZB for an amount of the order 
of J. The explicit S(uj) at T = was calculated in the 
'ladder approximation' by Elliott and Thorpe and found 
in very good agreement with experiments^. 

The finite temperature calculation of the two-magnon 
Raman scattering cross section in the ordered region, 
up to T^0.95Tn was performed by Balucani and 
Tognettitf, calculating the two-magnon propagator in 
the 'ladder approximation', taking also into account the 
damping and the temperature renormalization of the 
magnons at the boundary of the Brillouin zone 11 . The 
calculated spectra S(uj), at increasing temperatures, were 
found in very good agreement with the experimental 
ones£ and their characteristic parameters (peak and line- 
width) permitted to determine the temperature behav- 
ior of the frequency and damping of the ZB magnons 12 . 
In Figs, n and |21 we show the excellent agreement of 
our theoretical approach with the experiments in the or- 
dered phased. The validity of light scattering in prob- 
ing the characteristic of ZB magnons has been confirmed 
both from the theoretical and the experimental point of 
viewi^. In Fig. our theoretical ZB magnon damping 
calculations are compared with experimental data from 
different techniques. 

In the paramagnetic phase all experimental spectra 
show the persistence of a broad inelastic peak up to 
T ^\ATm ■ Only at T^Tm the spectra have a struc- 
tureless shape centered around uj = 0. As matter of fact, 
the highest wavevectors sample only the behavior of clus- 
ters of neighboring spins, thus giving a measure of the 
short-range antiferromagnetic order that is present at all 
finite temperatures. 

In the disordered phase conventional many-body meth- 
ods are of little use for a quantitative interpretation of the 
observed largely spread spectra. The concept of quasi- 
particle loses its meaning because of the overdamped 



FIG. 2: Experimental two-magnon spectra in KMF3 at dif- 
ferent temperatures-. 
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FIG. 3: Zone-boundary damping Tzb vs temperature. The 
symbols refer to different experimental approaches: in par- 
ticular the open circles are our light scattering data—. The 
dashed line is an improvement— to a previous (solid) theoret- 
ical curve—. 



character of the 'excitations'. The calculation of S(ui) 
can be instead approached by other more general the- 
oretical methods devoted to the representation of the 
dynamical correlation functions based on the linear re- 
sponse theoryii. Let us consider the 'Kubo relaxation 
function' associated with our scattering process: 

/o(<) - wmm Io dX { eAW ^(°)^ AW ^w) • («) 

Its Laplace transform fo(z) is related to the scattering 
cross section: 

g H°c i-e-fl" M/o(* = fo)- (9) 
Moriia has given the following continued fraction repre- 
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FIG. 4: Two-spin Stokes spectrum in KNiF 3 at T ~ 1.02 T N . 
The line reports the theoretical shape—, compared with ex- 
perimental data. 

sentation of the relaxation function^: 

/o(z)= z + A 1 / 1 (z) ; z + A n+1 /n+i(z)' 

( 10 ) 

which is formally exact, but allow us to do some approx- 
imations about the level of the termination f n+ i(z). The 
quantities A n can be expressed in terms of frequency mo- 
ments: 

/oo 
cW/oM. (ii) 
-oc 

In our calculations of fo(z) in the entire paramagnetic 
phas o 1 ^' 1 ^ , the coefficients Ai e A2 have been approx- 
imately evaluated by means of a decoupling procedure. 
Moreover, the third stage of the continued fraction 11 OB 
is evaluated assuming that 

A 3 / 3 (z)~A 3 [/3(0)+^(0)]. (12) 

The parameters involved in l|12f) can be estimated by the 
knowledge of the short time behavior of fo(t) determined 
by the first moments, (to 2 ) and (w 4 ). 

The results of our approach in the paramagnetic region 
are compared with the experiment in Fig. showing the 
persistence of the peak of the ZB magnetic excitations 
above the critical temperature. 



III. THE ONE-DIMENSIONAL 
ANTIFERROMAGNET 

A. Solitons in the antiferromagnet TMMC 

Interest in low-dimensional systems is motivated by 
the much greater simplicity of calculation as compared 
with the 3D ones. The powerful mathematical approach 
based on the inverse-scattering and Bethe Ansatz tech- 
niques permits to exactly solve some ID models, calcu- 
lating thermodynamic and sometimes transport quanti- 
ties both in classical and quantum casesit The most 



celebrated realizations of these models occur in ID mag- 
nets. An original suggestion by Mikeska 1 - was that the 
antiferromagnetic chain TMMC [(CH 3 ) 4 NMnCl 3 ] can be 
mapped onto a sine-Gordon classical ID model. The el- 
ementary excitations of the sine-Gordon field are given 
in terms of linear small-amplitude spin-waves and non- 
linear breathers and kink-solitons. The non- linear ele- 
mentary excitations give a detectable contribution to the 
magnetic specific heat. 

TMMC is composed of Heisenberg (S = 5/2) antiferro- 
magnetic chains along the z-axis: 

H = jJ2(SvS i+1 -6S?S* +1 ), (13) 

i 

with a very small easy-plane anisotropy (<5 = 0.0086). 

A magnetic field of the order of 1 10 T can be applied 
perpendicularly (y-axis) or along the chain . In the first 
case, with approximations the more valid the lower the 
magnetic field (H < 5 T) , in the continuum limit TMMC 
can be represented by the classical sine-Gordon Hamilto- 
nian: 

H=j J dx[i> 2 + cl<S>l + 2ul{l ~ cos *)] , (14) 

whose parameters are related with the magnetic Hamil- 
tonian H13[l. the reduced magnetic field h = g^H, and 
the lattice spacing a as follows: 



A= %Ta> c o = aJS ] /l--, u Q = hyjl--. (15) 

The energy of kink-soliton turns out to be 

E s = 8Auj c ~ h S , (16) 

and depends on the applied field. At difference with the 
ferromagnetic solitons, these solitons can be easily ex- 
cited at lowest temperatures and can give a significant 
contribution to the thermodynamics 19 . When the field 
is applied longitudinally along the z-axis only spin-waves 
are present: therefore, the specific-heat measurements 
were performed in the two configurations. The contribu- 
tion from the nonlinear excitations was obtained as the 
difference AC between the two experiments. 

The thermodynamic quantities were calculated by the 
classical transfer-matrix method^ for the sine-Gordon 
model <jl4"l) . We then used a classical discrete planar 
modelSi: 

H = J2 [2JS' 2 cos($ l -$ !+ i) + /iS*(l-cos$ l )] . (17) 

i 

verifying that it is qualitatively similar to the sine- 
Gordon. The comparison^ is shown in Fig. El where the 
linear spin-wave specific heat was subtracted to empha- 
size the nonlinear contribution, together with the predic- 
tion of the 'classical soliton gas phenomenology—. 
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FIG. 5: Experimental contribution of nonlinear excitations 
to the specific heat of TMMC, AC = C(H) - C(0) - ACsw- 
The field values are H = 5.39 T (•) and # = 2.5T (o). The 
dash-dotted line reports the result of the free soliton gas phe- 
nomenology. The planar model (interpolated crosses) appears 
to quantitatively explain the behavior of TMMC. 



This proved the presence of nonlinear excitations sim- 
ilar to sine-Gordon solitons, but the peak of the spe- 
cific heat occurs at temperatures where solitons can- 
not be considered to be non interacting and the 'classi- 
cal soliton-gas phenomenology' breaks down. When the 
magnetic field is increased up to 9.98 T the model l(T7|) 
is no more able to describe the experiments. A quasi- 
uniaxial model 21 was proposed and found in good agree- 
ment. For general reference on the subject see2£. 



B. The S—l quantum antiferromagnet 

We here deal with quantum antiferromagnetic spin 
chains, focusing our attention on the class of models de- 
fined by the Hamiltonian 



H 
J 



[(S*S* +1 +S?S? +1 )+\S?S? +1 +d(S?) 2 ] (18) 



with exchange integral J > and single-ion anisotropy d. 

One of the most surprising evidence of the differ- 
ence between ferro- and antiferromagnetic systems is re- 
lated with the so-called Haldane conjecture, i.e. with 
T = properties of integer-spin antiferromagnetic chains. 
In general, we expect three possible situations for the 
ground state of a magnetic system: either it is or- 
dered (with finitely constant correlation functions), or 
quasi-ordered (with power-law decaying correlation func- 
tions), or completely disordered (with exponentially de- 
caying correlation functions). One could intuitively ex- 
pect the third option to be possibly dismissed, based on 
the idea that, when thermal fluctuations are completely 
suppressed, the system be in an ordered or at least quasi- 
ordered ground state. This idea is in fact proved cor- 
rect for half-integer spin systems, thanks to the so called 



Lieb-Schultz-Mattis theorem 2 ^. Despite the generaliza- 
tion of such theorem to integer-spin systems being im- 
possible, its general validity has been taken for granted 
till 1983, when Haldanai suggested, for the integer- 
spin Heisenberg chain, an unexpected T = behavior: 
a unique and genuinely disordered ground state, mean- 
ing exponentially decaying correlation functions and a 
finite gap in the excitation spectrum. After more than 
two decades Haldane's idea that integer-spin systems can 
have a genuinely disordered ground state still stands as 
a conjecture. However, theoretical 2 ^* 2 ^ 2 ^!, experirnen- 
tal 28 1 29 1 30 1 3i 1 32 and nume rical 33 i 34 i 35 i 36 i 37 i 38 i 39 works have 
definitely confirmed its validity. 

Let us consider Eq. I|18l) for integer spin: in the (d, A) 
plane one may identify different quantum phases, corre- 
sponding to models whose ground states share common 
features. For A > three phases are singled out: the 
Neel phase (\^>d), where the ground state has a Neel- 
like structure, the so-called large-d phase (d^>X), where 
the ground state is characterized by a large majority of 
sites where S z — 0, and the Haldane phase, which extends 
around the isotropic point (<i = 0, A = l), and is charac- 
terized by disordered ground states. 

We first deal with the Ising limit, H/J = A £\ S?S? +1 : 
Upon its ground state, the antiferromagnetically ordered 
Neel state, one may construct three types of excitations: 
a single deviation, a direct soliton, an indirect soliton, 
where direct (indirect) refers to the fact that the excita- 
tion be generated by flipping all the spins on the right of 
a given site while keeping the z component of the spin 
on such site unchanged (setting it to zero). The above 
configurations have all energy +2A with respect to that 
of the ground state, and do generate, when properly com- 
bined, all the excited states; amongst them, we concen- 
trate upon those containing a couple of adjacent indirect 
solitons and notice that their energy is +3A, while excited 
states containing two separate indirect solitons have en- 
ergy +4A. Therefore, indirect solitons are characterized 
by a bounding energy A; moreover, one may easily see 
that isolated solitons may effectively introduce disorder 
in the global configuration of the system, while coupled 
solitons do only reduce the magnetization of each of the 
two antiferromagnetic sublattices^i. In fact, strings con- 
taining any odd (even) number of adjacent solitons act 
on the order of the global configuration as if they were 
isolated (coupled) solitons. 

As we move from the Ising limit, the transverse in- 
teraction J2i(Sf Sf+i + SfSf +1 ) comes into play, and is 
seen^i to more efficiently lower the energy of the system 
by delocalizing indirect solitons rather than single devi- 
ations or direct solitons, thus indicating configurations 
which uniquely contain indirect solitons as crucial in un- 
derstanding how the system evolves from the Ising limit 
(Neel phase) to the isotropic case(Haldane phase). 

From the above ideas we may draw a simple but sug- 
gestive scheme for such evolution: 

- in the Ising limit (A — > oo) the ground state is the anti- 
ferromagnetically ordered Neel state; 
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- as A decreases, indirect solitons appear along the chain 
in pairs, thus keeping the antiferromagnetic order; 

- as A is further lowered, indirect soliton pairs dissoci- 
ate due to the transverse interaction which, by spreading 
solitons along the chain, can cause the ground state to 
be disordered. 

Due to the privileged role of indirect solitons in the 
above scheme, we concentrate on configurations which do 
only contain indirect solitons. Such configurations gener- 
ate a subspace for the Hilbert space of the system, which 
is referred to, in the literature, as the reduced Hilbert 
spaced. States belonging to the reduced Hilbert space are 
strongly characterized by the fact that if one eliminates 
all sites with S z = 0, a perfectly antifcrromagnetically or- 
dered chain is left. Remarkably, this type of order, which 
is called hidden order in the literature, is not destroyed 
by soliton pairs dissociation, and it actually character- 
izes the disordered ground state of a Haldane system, as 
discussed below. 

In 1992 Kennedy and Tasaki (KT) defined a non-local 
unitary transformation^ which makes the hidden order 
visible, meanwhile clarifying its meaning. The transfor- 
mation is defined by U = (-l)^o+[iV/2] rj fe \j k w j t h 



i (V- £ t i s ? - i) ^ s t + 1 (> £j=i s ? +i 



where N is the number of sites of the chain, [iV/2] is 
the integer part of N/2, and N is the number of odd 
sites where S z = 0. If the pure state has hidden or- 
der, meaning that it only contains indirect solitons, then 
U\^) has spins with S z ^0 all parallel to each other. 
This point is made transparent by the introduction of 
the string order parameter— 



a 



string (ft) = . lim ( S i ex P 

|t-j|-KX5 



3-1 
ITT 53 5, 

l=i 



S?) , (19) 
H 



where a = x, y, z, and (•••)« indicates the expectation 
value over the ground state of the Hamiltonian TL. It may 
be shown that O z tling (TL) ^ if and only if the ground 
state belongs to the reduced Hilbert space. In other 
terms, while ferromagnetic order is revealed by the ferro- 
magnetic order parameter Of CTro = lim|j_ J -|_ >00 (5f S")-n , 
the hidden order is revealed by the string order param- 
eter Eq. (JTHJ. In fact, the non local transformation U 
relates the above order parameters through the relation 



oi^ = oi no (unu- 1 ) 



string 



(20) 



for r]=x,z, meaning that the analysis of the hidden order 
in a system described by TL may be developed by studying 
the ferromagnetic order in the system described by the 
transformed Hamiltonian TL=UTLU , which reads, for 
TL defined by Eq. |T5|). 



7 i=Y J ^S*S* +1 +S?e s *+ s ^S? +1 -\S z S z +1 +d(S z y 

(21) 



Our work developed as follows: one first assumes that 
the relevant configurations, as far as the Neel-Haldane 
transition is concerned, belong to the reduced Hilbert 
space; this permits, by the KT transformation, to re- 
strict the analysis to the subspace of states with either 
Sf = 1 or Sf =0, Mi. Then the expectation value of the 
transformed Hamiltonian TL Eq. l|21|l is minimized on a 
trial ground state whose structure takes into account 
at least short-range correlations between spins. By this 
procedure, we aim at following the effective dissociation 
of soliton pairs, in order to clarify the connection be- 
tween the occurrence of isolated solitons in the ground 
state, and the transition towards the completely disor- 
dered Haldane-phas o 39 i 41 ' 42 i 45 . 

In the framework of a standard variational approach, 
we should minimize ($o|ft|$o) with respect to a cer- 
tain number of variational parameters entering the ex- 
pression of the normalized trial ground-state |$o). By 
applying the non-local unitary transformation U we in- 
stead minimize {^ \UHU~ l *o) with |* ) = U\$ a ). 
and the transformed hamiltonian TL = UTLU^ 1 defined 
by Eq. J2U; if |$o) belongs to the reduced Hilbert space, 
it is 



|*o> =«|$o) 



y^c{s}isis 2 



sn) 



(22) 



with {s} = (s 1 ,s 2 ,s 3 ...s N ), and Si = (* |Sil*o} = +1, 0. 

The simplest trial ground state allowing the de- 
scription of soliton pairs dissociation is that defined 
by Eq. with c {s} = t SlS2S3 t S2S3Si ■ ■ ■ t SN _ 2SN _ lSN . 

The variational parameters are the six amplitudes 

t+++, *++o = *o++, i+o+: *o+o, *oo+=*+oo7 *ooo, 
where \t Si _ lSi s i+1 2 represents the probability for 
(Sf_i, S*, S z +1 ) to be equal (sj_i, Sj, Sj+i); a common ar- 
bitrary factor may be used for normalizing l^o)- We 
notice that the chosen form for c{ s } is such that the 
probability for |^ ) to contain coupled solitons is finite 
independently of that relative to the occurrence of iso- 
lated solitons, whose presence is unambiguously marked 
by t+o+ 7^0. 

Without going into the details of the variational cal- 
culations, reported in Ref. 14^. we here discuss our final 
results. Due to the normalization condition, the num- 
ber of variational parameters is reduced from six to five; 
moreover, the energy ( 1 >Sf \L{TLU~ 1 \ 1 $>o) is found to depend 
just on four precise combinations of the original param- 
eters, 



7 = |< ++0 | 2 |i o+| 2 

7T = l<++o| 2 |<+0+| 

X = I *oo+ 1 2 |*o+o I 

P = 1*000 1 



(••• + +00 + + •••) 

(••• + +0 + + •••) 

(• • • 00 + 00 • • • ) 

(• • • 000 • • ■ ) , (23) 



whose square moduli are related to the probabilities that 
the corresponding strings ( ) be contained in |\&o); m 
particular, 7 2 and it 2 refer to the probabilities for coupled 
and isolated solitons, respectively, to occur in the ground 
state. 
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FIG. 6: Phase diagram in the A > half-plane: our results 
(squares) are shown together with those of Ref. |43 (dotted 
lines); the Haldane phase should correspond to the shaded 
area, according to the best available numerical data£- (solid 
lines). 



Both the analytical expression for the energy and the 
numerical minimization show that it exists a critical value 
A c = A c (<i) > d such that, for A > A c the minimal energy is 
attained for 7r = p = 0; the condition A = A c (d) can hence 
define a curve of phase separation. We therefore single 
out three different phases, characterized by 

{a)ir = p = Q, (b) all parameters ^ 0, (c)x = ^+++ = 

in the ground state. The corresponding phase diagram 
is shown in Fig. together with that obtained with a 
factorized trial ground stated, and by numerical simula- 
tions^. 

The (a) — (b) transition is seen to quite precisely de- 
scribe the Neel— Haldane one, and this leads us to define 
the condition (a), meaning the occurrence of exclusively 
coupled solitons, as typical of the Neel phase. As for 
the (c) — (b) transition, it is to be noticed that the use 
of the reduced Hilbert space is not fully justified in the 
A < d region, where we in fact do not expect quantita- 
tively precise results. 

As for a comparison between our results and the exact 
numerical data available, we have considered, along the 
d = axis, two specific quantities: the critical anisotropy 
A c (d), where the Neel phase becomes unstable with re- 
spect to the Haldane one, and the ground-state energy 
Eo(d,\) at the isotropic point A=l. For the critical 
anisotropy we find A c (0) = 1.2044(5) to be compared with 
the value obtained by exact diagonalization^i, A c (0) s» 
1.19; for the energy we find E (0, 1) = - 1.3663(5) to be 
compared with Eq(0, 1) = — 1.4014(5), again from exact 
diagonalization technique^; the value obtained with the 
factorized trial ground state is^ Eq(0, 1) = — 4/3. 

In Fig. [7| we show the variational parameters as A is 
varied with e? = 0, i.e. along the y axis of the phase- 
diagram; in fact, rather than the parameters with respect 
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FIG. 7: Parameters (squares), W(2) (circles), wm (up- 
ward triangles), and u>(2,2) (downward triangles), for d = 0. 



to which we have actually minimized the energy, the fol- 
lowing combinations are considered: 

= tt 2 /4++ i w (2) = 7 ; 

t"(2,2) = (X 7 i+++) 1/2 5 W(3) = t 2 +++ ( P 7) 2/3 -(25) 

The above quantities have a straightforward physical 
meaning, as they are directly related with the probabili- 
ties for a soliton to appear along the chain as an isolated 
excitation (u>(i)), as part of a soliton pair (2w%\)i as part 
of a string made of three adjacent solitons (3w? 3 s), and 
finally as part of a string made of two soliton pairs sepa- 
rated by one site (4w* 2 2 j ) ■ From Fig. [3 it turns evident 
that the Haldane phase is featured by the occurrence of 
isolated solitons (wm ^0), as well as of strings made of 
three adjacent solitons (iu( 3 )^0). 

This result confirms that, as elicited by the analysis of 
the phase-diagram, the Haldane phase is characterized 
by our condition (6). 

Given their essential role, we have also studied the x 
and z component of the string order parameter, as well as 
the solitons density no = 1 — ((S* 2 ) 2 ). After KT we expect 
Og tring (W) ^ in both the Neel and the Haldane phase, 
and Og tring (7i) ^ just in the Haldane phase. In fact, 
analytical expressions for O x and O z may be written^ 
in terms of four of the five variational parameters 1)25(1 . 
and show that 

- Og tring (7i) = if 7r — p — or \ = t +++ = 0, i.e. in phase 
(a) and (c); 

- 0* tli (Ti) > in all phases, asymptotically vanishing 
as p— >1, i.e. in the far large-d phase. 

In more details, we notice that 0* tting — whenever 
the ground state does not contain strings made of an 
odd number of adjacent spins; as soon as the shortest 
string of such type, namely the isolated soliton, appears 
along the chain, then O x tring gets finite. The unphysi- 
cal result O z tTi > in the (c) phase, vanishing only as 
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FIG. 8: String order parameters Og tring (squares), OJ tr i ng 
(circles), and solitons density no (triangles), for d = 0. 



universality class the Haldane transition is suggested to 
belong to—. At the isotropic point (d = 0, A = 1) we find 
^string — 0.3700(5) in full agreement with the value ob- 
tained by exact diagonalization 48 . 

The overall good agreement between our results and 
the numerical available data, allows us to conclude that 
the Neel-Haldane transition is a second-order one, and 
that the string order parameter Og tring , revealing hidden 
order along the x direction, is the appropriate order pa- 
rameter for the Haldane phase. The disordered ground 
state featuring the Haldane phase is seen to originate by 
soliton pairs dissociation, according to this path: Solitons 
occur just in pairs in the antiferromagnetically ordered 
Neel phase; at the Neel-Haldane transition soliton pairs 
dissociate and the byproducts rearrange in strings made 
of an odd number of solitons. These strings are ulti- 
mately responsible for the disorder of the ground state. 
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FIG. 9: Critical behavior of 0^ tving for d = 0: squares 
are results; curves are obtained by best-fit procedure from 
Oatring ~ (A - \ c f with /3 fixed to 0.125 (dashed curve), and 
as fitting parameter, resulting in j3 = 0.217 (solid curve). Both 
procedures give A c = 1.2044(5), marked by a circle in figure. 



d — > oo rather than everywhere in the large-d phase, is 
due to our assuming the ground state to belong to the 
reduced Hilbcrt space, which is actually licit just in the 
A > d region. 

In Fig. |SJwe show Og tring , Of tring , and hq as A varies 
with d = 0: We underline that Og tring gets finite contin- 
uously but with discontinuous derivative at the transi- 
tion (reflecting the behavior of turn and wt^ shown in 
Fig-El, so that the Neel-Haldane quantum phase tran- 
sition is recognized as a second order one. In Fig. 
we zoom the order parameter Os tr j ng around the critical 
point: its behavior is seen to be described by a power law 
^string ~ (^c — ^)^> as expected for a continuous phase 
transition; our estimated value for the critical exponent is 
/3 = 0.217(5) to be compared with (3 = 0.125, correspond- 
ing to the Ising model in a transverse field, to whose 



IV. TWO-DIMENSIONAL ISOTROPIC 
HEISENBERG MODEL 

The 2D isotropic QHAF on the square lattice is one of 
the magnetic models most intensively investigated in the 
last two decades. This is due both to its theoretically 
challenging properties and to its being considered the 
best candidate for modeling the magnetic behavior of the 
parent compounds of some high-T c superconductors^. 

From a theoretical point of view the fully isotropic 
Heisenberg model in d dimensions, thanks to the simple 
structure of its Hamiltonian (whose high symmetry is re- 
sponsible for most of its peculiar features), may well be 
considered a cornerstone of the modern theory of critical 
phenomena, with its relevance extending well beyond the 
only magnetic systems. The d = 2 case earned additional 
interest, representing the boundary dimension separating 
systems with and without long-range order at finite tem- 
perature^. The antiferromagnetic coupling adds further 
appeal, as the classical-like Neel state is made unstable 
by quantum fluctuations and the ground state of the sys- 
tem is not exactly known. It can be rigorously proven^ 
to be ordered for 5>1; for 5 = 1/2 there is no rigorous 
proof, although evidences for an ordered ground state can 
be drawn from many different studies (for a review, see, 
for instance, Ref . l50]) . 

On the experimental side the attention on the proper- 
ties of 2D QHAF was mainly triggered by the fact that 
among the best experimental realizations of this model 
we find several parent compounds of high-T c supercon- 
ductors, as, e.g., La 2 Cu0 4 or S^CuC^C h, 51 ! 52 ! 53 , both 
having spin 5=1/2. In such materials, as well as in 
other magnetic compounds with a layered crystal struc- 
ture as La 2 Ni0^4 and KaNiF^S^ (5 = 1), Rb 2 MnF^ 
and KFeF^Si (5 = 5/2) or copper formate tetradeuter- 
ate (CFTD, 5=1/2)S the magnetic ions form parallel 
planes and interact strongly only if belonging to the same 
plane. The interplane interaction in these compounds 
is several orders of magnitude smaller than the intra- 
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plane one, thus offering a large temperature region where 
their magnetic behavior is indeed 2D down to those low 
temperatures where the weak interplane interaction be- 
comes relevant, driving the system towards a 3D ordered 
phase: an antiferromagnetic Heisenberg interaction and 
the small spin value make these compounds behave as 
2D QHAFs. Even the onset of 3D magnetic long-range 
order is however strongly affected by the 2D properties of 
the system: indeed, the observed 3D magnetic transition 
temperature is comparable with the intraplane interac- 
tion energy, i.e., several order of magnitude larger than 
that one can expect only on the basis of the value of 
the interplane coupling. Such apparently odd behavior 
can be easily understood by observing that the establish- 
ing of in-plane correlations on a characteristic distance £ 
effectively enhances the interplane coupling by a factor 
(£/a) 2 , a being the lattice constant. The latter consid- 
eration is one of the reasons explaining why most of the 
attention, both from the experimental and theoretical 
point of view, was devoted to the low-temperature be- 
havior of the correlation length £ of the 2D QHAF (in 
the following £ will be always given in units of the lattice 
constant a). 

The 2D QHAF is described by the Hamiltonian 

n=^J2Si-S i+d , (26) 

where J is positive and the quantum spin operators Si 
satisfy \Si\ 2 — S(S+1). The index i = (£1,12) runs over 
the sites of a square lattice, and d represents the displace- 
ments of the 4 nearest-neighbors of each site, (±1, 0) and 
(0,±1). 

In addition to the first approximations usually em- 
ployed to investigate the low temperature properties of 
magnetic systems as, e.g., mean- field and (modified) 
spin-wave theory, the critical behavior of the 2D QHAF 
was commonly interpreted on the basis of the results 
obtained by field theory starting from the so-called 2D 
quantum nonlinear a model (QNLctM) 58 , whose action is 
given by 

S=^-JdxJ^ dr (|Vn| 2 + |d T n| 2 ) ; \n\ 2 = 1 . (27) 

In the last Equation n(x) is a unitary 3D vector 
field, g = cA/p and u = cA/T are the coupling and the 
imaginary-time cut-off respectively, and the two param- 
eters p and c are usually referred to as spin stiffness and 
spin-wave velocity. Despite their names, the two param- 
eters p and c are however just phenomenological fitting 
constants which can be rigorously related to the proper 
parameter J and S of the original magnetic Hamilto- 
nian (|26l) only in the large-iS limiliiiS. The source of 
non-linearity in the model Eq. I|27|) . which is seemingly 
quadratic in the field variables, is the constraint imposed 
onto the length of the field n. 

The relation between the 2D QHAF and the QNLctM 
was first exploited to interpret the experimental data on 



cuprous oxides by Chakravarty, Halperin, and Nelson^ 
(CHN) who used symmetry arguments to show that the 
long- wavelength physics of the QHAF is the same of that 
of the QNLctM; in other words the physical observables 
of the two models show the same functional dependence 
upon T, if the long-wavelength excitations are assumed 
to be the only relevant ones, as one expects to be at low 
temperature. 

The analysis carried out by CHN on the QNLctM lead 
to single out three different regimes, called quantum dis- 
ordered, quantum critical (QCR) and renormalized clas- 
sical (RCR), the most striking difference amongst them 
being the temperature dependence of the spin correla- 
tions. If g is such as to guarantee LRO at T = 0, the 
QNLctM is in the RCR at very low-temperature and the 
correlation length £ behaves as§£: 

CHN found also that by raising the temperature any 
2D QNLctM with an ordered ground state crosses over 
from the RCR to the QCR, characterized by a correla- 
tion length £ oc a(T) = c/T. 

The first direct comparison between experimental data 
on spin 1/2 compounds and the prediction of the QNLctM 
field theory in the RCR gave surprisingly good agree- 
ment and caused an intense activity, both theoretical and 
experimental, in the subsequent years. However, with 
the accumulation of new experimental data on higher 
spin compounds it clearly emerged that the experimen- 
tally observed behavior of £(T) for larger spin could not 
be reproduced neither by the original simplified (2-loop) 
form of Eq. I|28|l given by CHN (which does not contain 
the term in square brackets), nor by the three- loops re- 
sult l|28|) derived by Hasenfratz and Niedermayer 60 (HN) ; 
moreover no trace of QCR behavior was found in pure 
compounds. The discrepancies observed could be due 
to the fact that the real compounds do not behave like 
2D QHAF or to an actual inadequacy of the theory. 
In particular the CHN-HN scheme introduces two pos- 
sible reasons for such inadequacy to occur: the physics 
of the 2D QHAF is not properly described by that of the 
2D QNLctM and/or the two(three)-loop renormalization- 
group expressions derived by CHN-HN do hold at tem- 
peratures lower than those experimentally accessible. Af- 
ter an almost ten years long debate, the latter possibility 
has finally emerged as the correct one, being strongly sup- 
ported not only by our own work, but also by other inde- 
pendent theoretical approaches^, joined with the anal- 
ysis of the experimental and the most recent quantum 
Monte Carlo (QMC) data for the 2D QHAF. 

The theoretical approach we employed to inves- 
tigate the 2D QHAF is the effective Hamiltonian 
methoo&SiSi, developed within the framework of the 
pure-quantum self- consistent harmonic approximation 
(PQSCHA) we introduced at the beginning of the 
gO's^ee. The PQSCHA starts from the Hamiltonian 
path-integral formulation of statistical mechanics which 
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allows one to separate in a natural way classical and 
quantum fluctuations: only the latter are then treated in 
a self-consistent harmonic approximation, finally getting 
an effective classical Hamiltonian, whose properties can 
thereafter be investigated by all the techniques available 
for classical systems. The idea of separating classical and 
quantum fluctuations turned out to be fruitful not only in 
view of the implementation of the PQSCHA, but also in 
the final understanding of the connection between semi- 
classical approaches and quantum field theories^, which 
could be possible also thanks to the paper by Hasen- 
frata§& about corrections to the field-theoretical results 
due to cutoff effects. 

The PQSCHA naturally applies to bosonic systems, 
whose Hamiltonian is written in terms of conjugate 
operators q = (q±, ...gjv), p = (pi,...pn) such that 
[Qm,Pn] =iS mn ; the method, however, does not require 
H(p,q) to be standard, i.e., with separate quadratic ki- 
netic p-dependent and potential ^-dependent terms, and 
its application may be extended also to magnetic sys- 
tems, according to the following scheme^: The spin 
Hamiltonian H(S_) is mapped to TL{p, q) by a suitable 
spin-boson transformation; once the corresponding Weyl 
symbol H(p,q), with p= (pi, ...p N ) and q = (q ± , ...q N ) 
classical phase-space variables, has been determined, the 
PQSCHA renormalizations may be evaluated and the ef- 
fective classical Hamiltonian Ti c s(j>, q) and effective clas- 
sical function O c ff(p,q) corresponding to the observable 
O of interest follow. Finally the effective functions 
7~l e s(jL) and O c fr(s), both depending on classical spin 
variables s with \s \ — 1 and containing temperature- and 
spin-dependent quantum renormalizcd parameters, are 
reconstructed by the inverse of the classical analogue of 
the spin-boson transformation used at the beginning. 

In order to successfully carry out such renormaliza- 
tion scheme, the Weyl symbol of the bosonic Hamilto- 
nian must be a well-behaved function in the whole phase 
space. Spin-boson transformations, on the other hand, 
can introduce singularities as a consequence of the topo- 
logical impossibility of a global mapping of a spherical 
phase space into a flat one. The choice of the trans- 
formation must then be such that the singularities oc- 
cur for configurations which are not thermodynamically 
relevant, and whose contribution may be hence approxi- 
mated. Most of the methods for studying magnetic sys- 
tems do in fact share this problem with the PQSCHA; 
what makes the difference is that by using the PQSCHA 
one separates the classical from the pure-quantum con- 
tribution to the thermal fluctuations, and the approxi- 
mation only regards the latter, being the former exactly 
taken into account when the effective Hamiltonian is re- 
cast in the form of a classical spin Hamiltonian. 

The spin-boson transformation which constitutes the 
first step of the magnetic PQSCHA is chosen according 
to the symmetry properties of the original Hamiltonian 
and of its ground state. In the case of the 2D QHAF 
both Dyson-Maleev and Holstein-Primakoff transforma- 



tion can be employed finally obtaining^: 

Jb i,d 

m = ^E 1 *^- 2 ^ (30) 

with the temperature and spin dependent parameters 

2 = l-j, (31) 

k 

fk = , C k = coth f k - i . (33) 
2 St Jk 

In the previous Equations 7& = (cosfci + cosfc2)/2, N is 
the number of sites of the lattice and k = (fci,^) is 
the wave- vector in the first Brillouin zone; S = S + 1/2 is 
the effective classical spin length, which naturally fol- 
lows from the renormalization scheme, and t = T/JS 2 
is the reduced temperature defined in terms of the en- 
ergy scale JS 2 . The renormalization scheme is closed by 
the self-consistent solution of the two coupled equations 
w fe =4K 2 (l - 7 2 ) 1 / 2 and k 2 = 6 2 - t/{2n 2 ). The pure- 
quantum renormalization coefficient T> = T>(S, t) takes the 
main contribution from the high-frequency part (short- 
wavelength) of the spin-wave spectrum, because of the 
appearance of the Langevin function £ k - V measures the 
strength of the pure-quantum fluctuations, whose contri- 
bution to the thermodynamics of the system is the only 
approximated one in the PQSCHA scheme. The theory 
is hence quantitatively meaningful as far as T> is small 
enough to justify the self-consistent harmonic treatment 
of the pure-quantum effects. In particular, the simple cri- 
terion T> < 0.5 is a reasonable one to assess the validity 
of the final results. 

The most relevant information we get from Eq. (|29|) is 
that the symmetry of the Hamiltonian is left unchanged 
so that from a macroscopic point of view the quantum 
system essentially behaves, at an actual temperature t, as 
its classical counterpart does at an effective temperature 
toff' = t/9 4 (S, t) . This allows us to deduce the behavior 
of many observables (but not all ! , see Refs. IbIbIbI for 
details) directly from the behavior of the corresponding 
classical quantities. This is the case of the correlation 
length, which turns out to be given simply by: 

£(*)=&l(teff) (34) 

so that once (9 4 (S, t) has been evaluated, the only addi- 
tional information we need is the classical £ c i(t), which 
is available from classical Monte Carlo simulation and 
analytical asymptotic expressions^, as t— »0. 

Sample results obtained by PQSCHA are shown in the 
figures. In Fig. ^]the correlation length for S—l/2 and 
S = 1 is compared with experimental data; a similar com- 
parison, including MC data for S = l/2 and experimen- 
tal data on S = 5/2 compounds KFeF4 and Rb2MnF4 is 
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FIG. 10: Correlation length £ vs t, for 5=1/2 (left- 
most) and 5 = 1. The symbols are experimental data; for 
5 = 1/2: 63 Cu NQR data£i (circles) and neutron scatter- 
ing data for La2Cu04 (squares—) and for Sr2Cu02Ci2 (up- 
triangles£2iS) ; for 5=1: neutron scattering data for La2Ni04 
(down-triangles—) and for K 2 NiF 4 (diamond a 52 ' 53 ). The clas- 
sical result (dash-dotted line) is also reported. 
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FIG. 11: The function y(t) —t ln£ vs t, for (from the right- 
most curve) 5 = oo, 5/2, 1 and 1/2; the up-triangles— and the 
diamonds^- are quantum MC data for 5 = 1/2; also reported 
are neutron scattering data for La2Ni04 (open circles—), for 
K 2 NiF 4 (squares^ 2 ^), KFeF 4 (filled circled) and Rb 2 MnF 4 
(down-triangles—). The abrupt rising of the experimental 
data for the 5 = 5/2 compounds at f~0.65 is due to the ef- 
fect of the small, but finite, anisotropies and will be discussed 
in more details in Section IVI 

made in Fig. El but along the vertical axis the quan- 
tity t In £ is reported in order to better appreciate the 
deviation from the predicted RCR behavior, that would 
correspond to a straight line at low t. From the last pic- 
ture one can easily see that both PQSCHA curves and 
experimental data for S > 1 (including the exact 5 = 00 
classical result) display a change of slope at intermediate 
temperature, followed by a curvature inversion at lower 
t. On the other hand, by looking at the S=l/2 case 
it becomes clear why the QNLctM approach gave such a 




FIG. 12: Correlation length £ and staggered susceptibility 
X* = x/S' 2 vs t for 5 = 1. Symbols are QMC data from Ref.[fl 



good agreement when firstly used to fit the experimental 
data. The change in both the slope and the curvature 
of t In £ is less pronounced and possibly occurs at lower 
temperatures, the lower the spin: in the 5=1/2 case, it 
is difficult to say whether these features are still present 
or not, but, if yes, they occur in a temperature region 
where the extremely high value of £ (~ 10 4 ) makes both 
the experimental and the simulation data more difficult 
to be obtained. 

After having realized that the field theoretical predic- 
tion by CHN could not be applied to the 5 > 1 2D QHAF 
in the temperature range probed by the experiments, the 
following questions were waiting for a satisfactory an- 
swer: (i) the real range of applicability of the asymp- 
totic three-loop expression (|28|l at different 5, and (ii) 
the possible extension of the PQSCHA results to lower 
temperature, both in view of (iii) a comprehensive de- 
scription of the behavior of the correlation length of the 
2D QHAF in the entire range of temperature and spin 
values. 

A substantial contribution to settle this conundrum 
came only from QMC simulations for higher spin values 
able to probe the very large correlation length region 73 : 
indeed, high precision Monte Carlo data for 5 = 1 and 
moderate correlation length could still be very well in- 
terpreted by PQSCHA and did not display the RCR 
asymptotic behavior, as shown^ in Fig. 1121 where we 
compared our curves for £ and staggered susceptibility 
X* with QMC data obtained by Harada et alw 

QMC results for £ by Beard and coworkers^ 3 - showed 
unambiguously that the three-loop Eq. (|28|l holds only 
for temperatures low enough to ensure an extremely large 
correlation length, e.g., £ > 10 5 for 5=1, £ > 10 12 for 
5 = 3/2, and generally cosmological correlation lengths 
for 5>3/2, thus definitely excluding any possibility of 
employing QNLoM results to interpret available experi- 
mental data. 

In Ref . l68l Hasenfratz showed why cutoff effects, which 
are so devious for 5 = 1/2, significantly modify the corre- 
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in its standard (low-T) version: the behavior of £ in the 
full temperature and spin- value ranges can thus be quan- 
titatively described by Eq. (|51Jl without any adjustable 
fitting parameter. 

The results obtained by the PQSCHA about the cor- 
relation length and other static quantities can also repre- 
sent the needed information to be inserted within other 
frameworks, like mode- coupling theory, to interpret ex- 
periments probing dynamic quantities, like nuclear mag- 
netic resonance (NMR): an example of the successful 
combination of PQSCHA and mode-coupling theory is 
given in Ref. 



V. TWO-DIMENSIONAL ANISOTROPIC 
HEISENBERG MODEL 



FIG. 13: Ratio f/f 3 l vs T/JS 2 for 5 = 5/2. Solid line: 
£d(£rff); dash-dotted line: £h — CeiI^St); dashed line: stan- 
dard PQSCHA result; symbols are QMC dataSi. 



lation length for S > 1 . Reaching such goal was possible 
only going back to a direct mapping between the QHAF 
and the QNLerM, and the resulting cutoff-corrected field- 
theoretical outcome is££ 



UT,S)=^(T,S) e 



-C(T,S) 



(35) 



where C(T,S), defined in Eq. (14) of Ref. M, is an in- 
tegral of familiar spin- wave quantities over the first Bril- 
louin zone. 

With this correction, which is the leading order in the 
spin- wave expansion for the cutoff correction, it is possi- 
ble to obtain numerically accurate agreement with QMC 
data down to £ > 10 3 for all S. 

In our most recent paper— about 2D QHAF we showed 
that by employing the explicit expression for C(T,S), 
and by substituting in Eq. I|28|l the spin stiffness p and 
the spin-wave velocity c given by the mapping of QHAF 
onto the QNLerM, the leading terms of the result (|35|l not 
only can be cast into the form £ H (T, S) — ^\(t^ s ), in strict 
analogy with the PQSCHA expression (|34(1 . but the effec- 
tive temperature is defined through a renormalization 
constant which is again a function of the pure-quantum 
fluctuations only! Such remarkable and unexpected fea- 
ture of the cutoff-corrected field-theory prediction, sug- 
gested us to substitute the perturbative expression £g} 
with the exactly known classical £ c i thus getting the re- 
sults presented in Fig. 1131 

It is thus made clear that the main features of the 
quantum correlation length at intermediate temperatures 
are due to essentially classical non-linear effects, which 
cannot be taken into account by perturbative approaches. 
Moreover, the effective exchange constant which defines 
the effective temperature tf s is seen to depend on the 
same pure-quantum renormalization coefficients defined 
by the PQSCHA, according to an expression which is 
very similar (equal) to that found by the latter approach 



While the theoretical debate mentioned in the previ- 
ous Section has been mainly dedicated to the isotropic 2D 
QHAF, real compounds are not actually well described 
by the isotropic model when the temperature is low: in- 
deed, the Mermin- Wagner theorem 4 states that a finite- 
temperature transition cannot occur in the 2D isotropic 
QHAF, while the experimental evidence of a transition 
suggests that 3D correlations and anisotropy effects, as 
well as a combination of both, must be considered. Easy- 
axis (EA) or easy-plane (EP) anisotropies turn out to be 
fundamental in the analysis of the critical behavior. 



A. 2D antiferromagnet with easy-axis anisotropy 

Several works (see Ref. for a review) have shown 
that many additional interaction mechanisms may be 
taken into account by inserting proper anisotropy terms 
in the magnetic Hamiltonian; in particular, the transi- 
tion observed in K 2 NiF 4 7e! {S = l), Rb 2 FeF^ (-5 = 2), 
K 2 MnF 4 77 , Rb 2 MnF4iS (5 = 5/2), and others, is seen 
to be possibly due to an easy-axis anisotropy. Such 
anisotropy has been often described in the literature 
through an external staggered magnetic field in order to 
allow for a qualitative description of the experimental 
data. However, this choice lacks the fundamental prop- 
erty of describing a genuine phase transition, as the field 
explicitly breaks the symmetry and makes the model or- 
dered at all temperatures. To preserve the symmetry 
under inversion along the easy-axis, it is actually appro- 
priate to insert an exchange anisotropy term in the spin 
Hamiltonian. Spontaneous symmetry breaking manifests 
then itself as a phase transition between ordered and dis- 
ordered states. The EA-QHAF Hamiltonian reads then 



H = 



S,; 



i+d 



qV qV \ i qz qz 



(36) 



where i = [ix, i 2 ) runs over the sites of a square lattice, d 
connects each site to its four nearest neighbors, J > is 
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the exchange integral and [i is the easy-axis anisotropy 
parameter (0 < fi < 1). Again, JS 2 = J(5+l/2) 2 sets the 
overall energy scale and t = T/JS 2 is the reduced tem- 
perature. When /i = 1 the model reduces to the isotropic 
QHAF. Note that a canonical transformation reversing 
the x and y spin components one one sublattice is equiv- 
alent to setting /i — > —jj,, so that the physical properties of 
the model are even functions of \i. The /z = case is called 
Ising limit, not to be confused with the genuine Ising 
model 78 , reproduced by Eq. (|36(l with fj, = and 5= 1/2. 
Despite being a very particular case of Eq. (|3^|) , the 2D 
Ising model on the square lattice is a fundamental point 
of reference for the study of the thermodynamic prop- 
erties of the EA-QHAF. A renormalization-group anal- 
ysis 7 ^ of the classical model predicted the occurrence of 
an Ising- like transition at a finite temperature 1%r(li) of 
the order of unity for any value of fi, no matter how near 
to the isotropic value jii = l; this analysis received the 
support of several Monte Carlo simulations£2i2ii§2i2ii. 

As for the quantum case, up to a few years ago no 
information was available about the value of the criti- 
cal temperature t c (/i, 5) as a function of anisotropy and 
spin, save the fact that t c (0, 1/2) =0.567J (Onsager so- 
lutionS 8 .) and t c (l, 5) =0 (isotropic limit). As a conse- 
quence it was also uncertain whether or not the small 
anisotropy (1 — /i~10~ 2 ) observed in real compounds 
could be responsible of transitions occurring at critical 
temperatures of the order of J, also accounting for the 
fact that quantum fluctuations are expected to lower the 
critical temperature with respect to the classical case. 

Over the last few years our group developed a quantita- 
tive analysis of several thermodynamic properties of the 
model, by means of the effective Hamiltonian methoo&S& 
for spin 5 > 1 and by means of quantum Monte Carlo 
simulations 84 in the case 5=1/2. 

The effective Hamiltonian 85 i 86 i 87 for the EA-QHAF is 
expressed for classical spins as 

= '-f E + S y A +d ) + 4s! +d ] , (37) 

Jt} i,d 

and shows a weaker renormalized exchange j c s(t,fi) < 1 
and easy-axis anisotropy fj, e g(t,fi) > n, besides an addi- 
tional free-energy term that is not reported. By means 
of TC c s a series of thermodynamic quantities were stud- 
ied2£i2&: internal energy, specific heat, staggered mag- 
netization, staggered correlation function, staggered cor- 
relation length, staggered susceptibility. This required 
extensive classical Monte Carlo simulations, as varying 
the temperature gives an effective system with different 
effective anisotropy /i e ff(i,/i). The quantum phase di- 
agram reported in Fig. ^] could be built upSLS in a 
simpler way starting from the knowledge of the classical 
one and using a relation that follows from the form of 
Eq. 123, 

t c (n, 5) = j eS (t, n, 5) t?(v eS (t, fi, 5)) . (38) 
In the region of very weak anisotropy, which is the most 
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FIG. 14: Critical temperature vs anisotropy fi for the 2D 
easy-axis antiferromagnet. Dotted line: the fit of ^(m) build 
up from classical MC data—*— *— (open triangles). Solid 
lines: PQSCHA result for 5=1 and 5/2. Quantum MC 
data (open circles— and diamonds—) for S — 1/2, asymptot- 
ically described by t c (fJ.) =2.49/ ln[70/(l-j«)] (dashed line), 
and exact result for the Ising model, fj, = (open square) . 
Experimental data for the compounds YBa2Cua06.1 (down 
triangleSi), K 2 NiF 4 (cross22^), Rb 2 Ni F 4 (up triangle^), 
Rb 2 MnCl 4 (circle 95 ), Rb 2 MnF 4 (diamond^ffi). In the inset 
the region of weak anisotropy is enlarged. 



important in view of the characterization of experimen- 
tally accessible materials, we verified that the Ising-like 
transition temperature decreases very slowly (logarithmi- 
cally) towards its vanishing value in the isotropic limit, 
so that t c remains substantially of the order of unity. 

As a sample of the various results that were obtained, 
we report in Fig. 1151 the comparison 8 ^ with the experi- 
mental dataffi for the correlation length of the 5 = 5/2 
magnet Rb2MnF 4 , that results quite well described by 
the anisotropic model with J = 7.42 K and /i = 0.9942. 
Rb2MnF 4 is known to behave as a 2D magnet both above 
and below the observed transition 7 ^, so that the critical 
behavior is not contaminated by the onset of 3D order 
and a clean characterization of the transition is possi- 
ble. In Ref. we have compared our theoretical re- 
sults also with the neutron scattering experimental data 
for the staggered magnetization, staggered susceptibility 
and correlation length of Rb2MnF 4 and found an excel- 
lent agreement both for the overall temperature behavior 
and for the value of the critical temperature, that per- 
fectly coincides with the one deriving from the experi- 
mental analysis, T c = 38.4 K (i.e., i c = 0.575). 

Another quantity that shows a signature of the an- 
isotropy and of the Ising transition is the specific heat: 
in Fig. 1161 a comparison with experimental data is shown 
in the case of the 5 = 5/2 compound Mn- formate di- 
Urea™, whose anisotropy can be estimated from the sole 
knowledge of the exchange integral and of the measured 
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FIG. 15: Correlation length vs t for S = 5/2, n = 0.9942 (full 
curve) and /i = l (isotropic, dash-dotted curve); the symbols 
are neutron scattering data- 7 for Rb2MnF4. The triangles are 
quantum Monte Carlo data— for the isotropic model. 
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FIG. 17: Uniform susceptibility of the EA model for fi — 0.99 
from quantum Monte Carlo simulations 84 . Full diamonds: 
longitudinal branch; open diamonds: transverse branch; stars: 
data for the isotropic model—-. The lines are guides for the 
eye. The arrow indicates the estimated critical temperature. 
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FIG. 16: Specific heat vs t = T/JS 2 , for 5 = 5/2. Mn- 
f-2U experiments (squares)iffi, EA-QHAF with ^ = 0.9942 
(solid line)^ 8 -, isotropic QHAF (dashed line). Note that 
the correct anisotropy for this compound is estimated to be 
fj, = 0.981. Error bars are due to the experimental uncertainty 
on J for Mn-f-2U. 



transition temperature to be /i = 0.981. The compar- 
ison reveals the existence of a crossover from a high- 
temperature 2D-Heisenberg regime to a critical 2D-Ising 
regime that triggers the observed 9 ^ 3D phase transition 
atT N =3.77K. 

Finally, for the strongest quantum case, 5=1/2, we 
have used the continuous-time quantum Monte Carlo 
method based on the loop algorithm— The general out- 
come of the numerical simulations is that the thermo- 
dynamics of 2D quantum antiferromagnets is extremely 
sensitive to the presence of weak easy-axis anisotropies 
of the order of those of real compounds. For instance, 
in Fig. El it is shown that for /i = 0.99 the uniform sus- 



ceptibility, which is a noncritical quantity, undoubtedly 
shows a characteristic anisotropic behavior with a differ- 
ent temperature dependence of the transverse and lon- 
gitudinal branches: the former displays a minimum and 
the latter monotonically goes to zero, as expected for 
an EA antiferromagnet. This behavior results from the 
anisotropy-induced spin ordering, that makes the system 
more sensitive to the application of a transverse mag- 
netic field, rather than of a longitudinal one. Both the 
minimum of the in-plane component and the decrease of 
the longitudinal one are close to the transition, a feature 
also peculiar to the Ising model. Results for the critical 
temperature at 5 = 1/2 are already included in Fig.lT4l 

B. 2D antiferromagnet with easy-plane anisotropy 

In the case of an easy-plane anisotropy the Mermin- 
Wagner theorem still holds, so that no finite-temperature 
transition towards a phase with a finite order parameter 
may occur. However, a Berezinskii-Kosterlitz-Thouless 
(BKT) transition 102 ! 103 , related with the existence of 
vortex-like topological excitations, is known to charac- 
terize the class of the easy-plane models and may occur 
at a critical temperature t c (X, 5) > 0. The reference sys- 
tem for the easy-plane class is the planar rotator model, 
or XY model, defined in terms of two-component spins: 
above t c the system is disordered, with exponentially de- 
caying correlation functions; in the region < t < t c the 
system is in a critical phase with vanishing magnetization 
and power-law decaying correlators (quasi- long-range or- 
der); at t — the magnetization gets a finite value and 
the system is ordered. 

The observation of clear signatures of BKT critical 
behavior in real magnets is a controversial issue. How- 
ever, it can explain the properties of several layered com- 
poundsSi 1 i24i 1 i2£* 1 i ii whose high-temperature phase can be 
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described by a purely 2D Heisenberg Hamiltonian, with 
an exchange interaction often displaying weak easy-plane 
(EP) anisotropies, on the order of 10~ 2 ~ 10~ 4 times 
the dominant isotropic coupling 105-106 . Symmetry and 
universality arguments suggest that the EP anisotropy 
drives the system towards a BKT behavior at finite tem- 
perature, and the enhanced intraplane correlations trig- 
ger the transition to the observed 3D ordered state. As 
a consequence, 2D critical behavior in close proximity 
of the would-be BKT transition is masked by these 3D 
effects. 

The EP-QHAF Hamiltonian reads 



H 



-E 

2 f-f 



S: 



s?s? +d + xs?s„ 



i+d 



(39) 



where A is the easy-plane anisotropy parameter 
(0 < A < 1). Again, JS 2 = J(5+l/2) 2 sets the overall en- 
ergy scale and t~T/JS 2 is the reduced temperature. 
When A = 1 the model reduces to the isotropic QHAF. 
Note that a canonical transformation reversing the x 
and y spin components one one sublattice is equivalent 
to setting (J, A) —* (— J, — A), so that negative values of 
A (— 1 < A < 0) correspond to the easy-plane ferromag- 
net. The A = case is called XY model or XXO model. 
However, at variance with the planar rotator model, 
out-of-plane fluctuations are present both in the clas- 
sical and in the quantum EP models. Nevertheless, if 
|A| < 1, the classical EP model still undergoes a BKT 
phase transition— £. MC simulations of the classical sys- 
tema 108 i 109 i 110 i ul confirm that the planar and the XXZ 
model share the same qualitative behavior, but the value 
of the transition temperature of the planar models low- 
ers by 22% in the XXO (i.e., with A = 0) model 110 i m , 
as a consequence of out-of-plane fluctuations. A renor- 
malization group calculation 107 predicts that the transi- 
tion temperature vanishes logarithmically as the isotropic 
limit |A| — > 1 is approached, and this was also verified in 
classical MC simulations 110 . Experiments 5 ^* 5 ^ and quan- 
tum MC simulations 113 indicated that the qualitative be- 
havior of the BKT transition is preserved in the quantum 
system, with only quantitative modifications of the crit- 
ical parameters due to the quantum fluctuations. 

65.66 



We applied the effective Hamiltonian formalism— *2ii to 
the EP-QHAF, finding that it was necessarj^ 7,114 to re- 
sort to the Villain or to the Holstein-Primakoff trans- 
formation, depending on the anisotropy being strong or 
weak, respectively. While the above approach gives reli- 
able results (with a smooth enough connection at inter- 
mediate anisotropy) for spin S > 1, we adopted quan- 
tum Monte Carlo simulationa 84 i 88 i 115 i 116 i 117 in the case 
£=1/2. The effective Hamiltonian 87 i 114 i 118 for the EP- 
QHAF, in terms of classical spins, takes the form 



JS 2 



(40) 



and displays a weaker renormalized exchange 
j e s{t,X)<J and easy-plane anisotropy X c s(t, A) > A 
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FIG. 18: BKT critical temperature vs anisotropy A for the 
easy-plane model with 5 = 1/2, 1, 5/2, oo. The triangles re- 
port the classical (5 = oo) MC data—- used to construct the 
curve for ^'(A) and hence the renormalized curves for 5 = 1 
and 5/2: the curves obtained through the Holstein-Primakoff 
(A<0.5) and the Villain (A > 0.5) spin-boson transforma- 
tion are seen to connect in a fairly smooth way. The di- 
amonds are our QMC data— for 5 = 1/2, while the circles 
are earlier QMC results—-. The inset expands the nearly 
isotropic region, in which the expected logarithmic behavior 
t c (A) ~ (a — In |1— A|) _1 is fitted by the dashed curve. 



(an additional free-energy term is not reported). In 
analogy to the EA case, the BKT transition temperature 
can be obtained by renormalization of the classical one 
using the self-consistent relation 

t c (X, S) = j eS {t, A, S) t?(\ eS (t, A, S)) . (41) 

In Fig. |^| the phase diagram of the EP-QHAF is re- 
ported, including the QMC results for 5 = 1/2. It is seen 
that the BKT transition temperature stays large (i.e., 
comparable to the exchange constant) also for very weak 
EP anisotropy. 

However, as explained above, the problem of detecting 
the incipient BKT transition requires to look for signa- 
tures of XY behavior in a region above the transition. We 
have showniii 5 " that a suitable quantity is the uniform sus- 
ceptibility x™, which has in-plane (a = x, y) and out-of- 
plane (pt = z) components and is noncritical, i.e., it does 
not show singularities at t c . Fig.lT§lshows indeed that %u 2 
deviates from the isotropic Xu and displays a minimum. 
A similar feature is also present in other quantities 84 ' 115 
and occurs around the temperature t co (X) that can be 
generally defined as the minimum of Xu Z (t, A). The pro- 
nounced deviation of Xu* from the isotropic behavior is 
due to a simple statistical reason. The uniform suscepti- 
bility arises from spin canting: two antiferromagnetically 
coupled spins in an infinitesimal magnetic field h mini- 
mize their energy when they lie almost orthogonal to h 
and slightly cant in its direction, thus giving a linear re- 
sponse; if they are locally parallel to h the response is 
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FIG. 19: Out-of-plane (zz) and in-plane (xx) uniform sus- 
ceptibility Xu for A = 0.999. Diamonds are our QMC results, 
crosses are experimental data4A2*12 £ for Sr2Cu02Cl2- The cir- 
cles report the result for the isotropic QHAF. The vertical 
lines mark the 3D transition temperature t N =0.176 and the 
crossover temperature t co — 0.227. 



instead negligible. When for t<t co the anisotropy be- 
comes effective, the fraction of spins aligned in the EP 
rapidly increases compared to that of the isotropic case 
(~ 2/3) and the response to a field along z (i.e., Xu Z ) i s 
proportionally larger. 

The layered cuprate Sr2Cu02Cl2 is a good realiza- 
tionS of the EP-QHAF model, with J=1450K; con- 
sidering the spin-wave gap renormalization, its bare 
anisotropy is estimated to be A ~ 0.0014. Experimen- 
tal data for the uniform susceptibility of this com- 
poundii2i±211) are reported in Fig. ED They excellently 
compare with our results at A = 0.999: the position of the 
minimum of \u Z gives t co — 0.227(15). Close to the crit- 
ical region the experimental data are affected by the 3D 
nature of the ordering of the real magnet: the Neel tran- 
sition is observed at t N =0.176(10) and well compares 
with the theoretical estimate t c (X=0. 0014) = 0.179(10), 
confirming that 3D ordering is induced by the incipient 
intra-layer BKT transition. 

It is worthwhile to mention that on the triangular lat- 
tice the easy-plane antiferromagnet has very peculiar be- 
havior, already at the classical level, due to the frus- 
tration effect of accommodating three antifcrromagnetic 
spins on a plaquette. Indeed the minimum energy cor- 
responds to a configuration with the three sublattices 
aligned in the easy plane at equal angles 27r/3. As clock- 
wise and counterclockwise plaquette vorticities are pos- 
sible, this configuration is twofold degenerate and there 
is chiral symmetry, which corresponds to an Ising-like 
order parameter. Therefore both a BKT and an Ising 
transition coexist in the system. We have studied the 
triangular antiferromagnet both in the classical^Si and 
in the quantum?^ case, constructing the phase diagram 
for varying anisotropy and showing that the transitions 
occur at slightly different temperatures. 



C. 2D antiferromagnet in an applied Zeeman field 

An interesting behavior is shown by the 2D Heisenberg 
antiferromagnet when a magnetic field is applied. Indeed, 
a frustration phenomenon occurs, as antiferromagnetism 
tends to antialign spins while the field tends to align them 
with itself: in the classical system^ it appears that the 
minimum-energy configuration is the one with the spins 
almost orthogonal to the field and canted in its direc- 
tion. Therefore, provided the field is not strong enough 
to overcome the exchange and to saturate the magneti- 
zation, it acts as an effective easy-plane anisotropy and 
one expects to observe BKT behavior. Remarkably, as 
this can be induced in a real quasi-2D antiferromagnetic 
system by means of an applied field, the strength of the 
effective anisotropy is in this case tunable. Even though 
2D criticality just acts as a trigger for 3D ordering, by 
observing that the critical temperature follows the pre- 
dicted behavior upon the field, an experiment could re- 
alize an objective observation of genuine 2D behavior. 

The 2D QHAF in a uniform magnetic field is described 
by the Hamiltonian 



9(J>, 



(42) 



where H is the applied Zeeman field, ^ B the Bohr mag- 
neton, and g the gyromagnetic ratio. 

We have studied 124 ! 125 the S = ± 2D QHAF in an uni- 
form magnetic field by means of the QMC method based 
on the worm algorithm. Our results confirmed that an 
arbitrarily small field is able to induce a BKT transition 
and an extended XY phase above it, as in the case of an 
easy-plane exchange anisotropy. The field-induced XY 
behavior becomes more and more marked for increasing 
fields, while for strong fields the antiferromagnetic behav- 
ior along the field axis is nearly washed out, so that the 
system behaves as a planar rotator model with antifer- 
romagnetism surviving in the orthogonal plane only; the 
BKT critical temperature, as reported in Fig 1251 (where 
t = T/J and h = 2g[i g H/ J), vanishes as the field reaches 
the saturation value h c and the effective rotator length 
goes to zero. We have therefore shown that the model in 
a moderately strong field represents an ideal realization 
of the XY model and that XY behavior can be detected 
by measuring standard non-critical quantities, as the spe- 
cific heat or the induced magnetization; an experimental 
realization of the XY model in purely magnetic systems 
and a systematic investigation of the dynamics of vor- 
tex/antivortex excitations is therefore possible. 



VI. FRUSTRATION IN THE 2D QUANTUM 
J1-J2 HEISENBERG MODEL 

The study of frustrated quantum spin systems is one 
of the most challenging and exciting topics in theoreti- 
cal magnetism because of the possible existence of a non 
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FIG. 20: Phase diagram of the S = § 2D QHAF in a magnetic 
field. Open symbols refer to the classical limit of the model—-; 
the triangle^ and the squares^ are QMC results. 



magnetic zero temperature phase. A very extensively 
investigated, yet largely debated model is the so-called 
J1-J2 Heisenberg model with competing antiferromag- 
nctic couplings (Ji,J2>0) between nearest-neighbors 
(nn) and next-nearest-neighbors (nnn) 



H = J\ Si-Sj + J2 Si-Sj 



(43) 



where the spin operators are defined on a periodic lat- 
tice with N = Lx L sites; hereafter a — J-ij J\ defines the 
frustration ratio. 

In the classical limit (S—>oo), the minimum energy 
configuration has conventional Neel order with magnetic 
wave vector Q = (it, it) for a < 0.5. Instead, for a > 0.5, 
the antiferromagnetic order is established independently 
on the two sublattices, with the two staggered magne- 
tizations free to rotate with respect to each other. One 
of the two families of collinear states, with pitch vectors 
Q = (ir,0) or (0,7r), are selected by an order-by-disorder 
mechanism as soon as thermal or quantum fluctuations 
are taken into account. As a result, for a > 0.5 the classi- 
cal ground state breaks not only the spin rotational and 
translational invariance of the Hamiltonian - as the con- 
ventional Neel phase - but also its invariance under tt/2 
lattice rotations, the resulting degeneracy corresponding 
to the group 0(3) x Z-i. Remarkably, the additional dis- 
crete Zi symmetry can in principle be broken at finite 
temperatures without violating the Mermin- Wagner the- 
orem. On this basis, in a seminal paper 1 - 7 , Chandra, 
Coleman, and Larkin (CCL) proposed that the 2D J1-J2 
model could sustain an Ising phase transition at finite- 
temperature, with an order parameter directly related to 
the Zi degree of freedom induced by frustration. They 
also provided quantitative estimates of the critical tem- 
peratures in the large-a limit for both the classical and 
the quantum cases. 

This transition in the classical model has been estab- 
lished by an extensive Monte CarloiSS. In the quantum 



case, the occurrence of a low-temperature phase with a 
discrete broken symmetry has been subject of debates 
in connection with the discovery of three vanadate com- 
pounds (Li 2 VOSi0 4 , Li 2 VOGe0 4 , and VOM0O4) whose 
relevant magnetic interactions involve nearest and next- 
nearest spin-1/2 V A+ ions on weakly coupled stacked 
planes. In particular, NMR and /xSR measurements on 
Li2VOSi04 129 indicate the occurrence of a transition to a 
low-temperature phase with collinear order at T N ~ 2.8 K. 
However in the experiments with vanadate compounds, 
structural distortions, interlayer and anisotropy effects 
are likely to come into plaj*i2!i, and on the other hand 
the theoretical investigation cannot rely on the insight 
provided by quantum Monte Carlo methods as their re- 
liability in presence of fru stration i s strongly limited (see 
the review articles in Ref. llSlllsi . 

A complete study of the thermodynamic properties of 
the the quantum J1-J2 model in its collinear phase has 
been pursued within the PQSCHA scheme^ described in 
Section IIVI by which the thermodynamics is rephrased 
in terms of a classical effective Hamiltonian with renor- 
malized parameters depending on the spin value, tem- 
perature, and frustration. It is possible to show that, to 
0(1/5), the effective Hamiltonian can be recast in a form 
preserving all the symmetries of the original model, and 
that reads (except for uniform terms): 



jfS 2 ^ 



(44) 



where s j are classical vectors of length 1 , S = S + i 
is the effective spin length, and jf = ((%+ B 2 ) Q\ Ji/2, 
J^ ff = 0|j 2 , are the quantum-renormalized exchange 
integrals, with spin-, temperature- and frustration- 
dependent renormalization parameters 9 X , 9 y and 62- 

The occurrence of the transitioniSJ. in the quantum 
case can be directly addressed within our approach by 
calculating the critical temperatures as functions of the 
spin and of the frustration ratioAS. Using a simple scal- 
ing argument the critical temperatures in the quantum 
case T C (S, a) can be related to those of the classical 
model Tc Cl \a) through the following self-consistent re- 
lation 87 ! 122 

T C (S, a) = jf(T c , S, a) T^(a cS (T C: S, a)) , (45) 

where jf = jfS 2 /Ji and a cS = Jf/Jf. The classi- 
cal transition temperature, Tc Cl \a) is accurately known 
through extensive MC simulations for a < 2; it vanishes 
for a — > 1/2 and grows more or less linearly for a > 1. 

The behavior of the transition temperature versus the 
frustration ration is plotted in Fig. for different val- 
ues of the spin length. In order to represent the whole 
interval of a £ [1/2, 00) in Fig. we have plotted both 
the MC and the CCL estimates of the classical critical 
temperatures as a function of a/(l + a). The mismatch 
between the MC and CCL predictions is a minor flaw 
that can be corrected by slightly modifying CCL's crite- 
rion for the determination of the transition temperature 
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FIG. 21: Renormalized critical temperature of the CCL tran- 
sition for vario us va lues of the spin. Classical data (•) are 
taken from Ref. Il28t The solid lines on the right are the CCL 
and the CCL* prediction for the classical and 5 = 1/2 case 
(see text). The ellipse on the horizontal axis mark the non- 
magnetic (spin-liquid) phase between a c (S — oo) = 0.5 and 
Oo(5=l/2)~0.6. 

as explained in Ref. Il28l Remarkably, while for large 
a the transition temperature vanishes for a — ■+ oo for any 
value of the spin, in the opposite limit the critical temper- 
atures vanish approaching a critical value a c > 0.5 that 



increases as S decreases, thus confirming the existence of 
a non-magnetic phase in the regime of high frustration. 
In particular for S= 1/2, a c ~0.6 in agreement with the 
previous estimates of the zero-temperature quantum crit- 
ical poin1ii2£. 



VII. CONCLUSIONS 

The activity in magnetism of the Condensed Matter 
Theory group in FirenzeA24 stems from the early work 
on two-magnon Raman scattering in the seventies, and 
has grown up in the years with the collaboration of sev- 
eral scientists. In this paper we summarize the relevant 
theoretical work that concerns antiferromagnetic mod- 
els. This activity has been mainly concentrated on low- 
dimensional systems and has found one of its main mo- 
tivations in the intent of interpreting the data collected 
in experiments on real materials. Among the prominent 
subjects, we reported about soliton-excitation effects in 
one-dimensional systems, critical and near-to-critical be- 
haviors and phase transitions in two-dimensional ones. 
Besides that, we also faced some intriguing theoretical 
problems where fundamental aspects of Quantum Me- 
chanics come into play, as, e.g., the ground state of 
antiferromagnetic chains with integer spin or the pos- 
sible quantum critical regime predicted from the field- 
theory treatment of the two-dimensional antiferromag- 
netic Heisenberg model. 
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